Theory of the field-induced gap in S = 1/2 antiferromagnetic chains 



Masaki Oshikawa 1 and Ian Affleck 1 ' 2 
Department of Physics and Astronomy 1 and Canadian Institute for Advanced Research 2 , 
University of British Columbia, Vancouver, BC, V6T 1Z1, CANADA 

(22 May 1997) 



OS 

o> 

=3 



i 

-4-* 



I 

C 

o 
o 



> 
in 

00 

o 

so 
O 

Os 

c3 



O 
o 



X 



In a recent neutron-scattering experiment on the quasi- 
one-dimensional 5=1/2 antiferromagnet Cu Benzoate, a 
gap was induced by an applied magnetic field. We argue 
that the primary mechanism of the gap formation is an ef- 
fective staggered field due to both the alternating g-tensor 
and the Dzyaloshinskii-Moriya interaction. We explain the 
dependence of the gap on the applied field, as well as identify 
several peaks in the structure factor S(q,u>). 

PACS number: 75.10Jm 



Quantum spin chains have attracted much interest for 
a long time. This is partly because sophisticated theo- 
retical analysis, such as exact solutions, can be applied 
to one dimensional systems. Not only are they easier 
to analyze, it has been also recognized that the effect of 
quantum fluctuations is more significant than in higher- 
dimensional systems, resulting in many interesting phe- 
nomena. On the other hand, progress in experimen- 
tal techniques has increased the opportunity to observe 
physics of one-dimensional systems. In a recent high field 
neutron-scattering experiment jjj on Cu Benzoate, which 
is a (quasi-)one-dimensional 5 = 1/2 antiferromagnet, 
the field-induced shift in the soft-mode momentum is ob- 
served for the first time. Although the shift of the mo- 
mentum is consistent with previous theoretical analysis 
on the Heisenberg antiferromagnetic chain, the experi- 
ment also found an unexpected excitation gap induced 
by the applied field. The observed gap is proportional to 
Hq ' 65 where Ho is the magnitude of the applied field. 
While the data is consistent with the power law with 
the same exponent 0.65 for three different directions of 
the applied field, the coefficient depends on the direction. 
The ratio of the coefficient is found to be 0.55 : 1.0 : 2.0 
for the field applied in a", b, c" axes, which are the prin- 
cipal axes of the effective exchange interaction. (For de- 
tailed description of the compound, see Ref. S.) The 
observed gap can be as large as 0.3 J where J is the ex- 
change coupling in the chain direction, at Hq = 7T where 
the average magnetization is 0.06 per site. 

In this letter, we discuss the mechanism of the field- 
induced gap observed in Cu Benzoate. We argue that the 
primary mechanism is due to an effective staggered field. 
As pointed out by Dender et al. , an effective staggered 
field is generated by the alternating g-tensor. We found 
that the effective staggered field is also generated by the 
Dzyaloshinskii-Moriya (DM) interaction, and the latter 



is no less important than the former. Our theory suc- 
cessfully explains the experimental data, including the 
angle dependence of the gap. For quasi-one-dimensional 
compounds with alternating crystal axes, both effects are 
expected. Thus, our theory should apply to such com- 
pounds in general. 

Since the interchain coupling is very weak in the com- 
pound, the gap formation should be understood primar- 
ily in a one-dimensional model. In this letter, we restrict 
our discussion to one-dimensional models. As a first ap- 
proximation, the system would be described as the stan- 
dard isotropic Heisenberg antiferromagnet. Actually, the 
neutron scattering data at zero magnetic field is consis- 
tent with the theoretical analysis based on the standard 
Heisenberg model with exchange coupling J = 1.57meV. 
However, the standard Heisenberg model in an applied 
field remains gapless from zero magnetic field up to satu- 
rated magnetization Thus we have to consider some 
modification. 

Even if we generalize the model Hamiltonian, the sys- 
tem remains gapless for generic values of the magnetiza- 
tion, as long as it has rotation symmetry about the di- 
rection of the magnetization (we call this axial symmetry 
hereafter). This can be seen from abelian bosonization 
or the generalized Lieb-Schultz-Mattis theorem j|. Since 
the gap is observed with a continuously changing mag- 
netization, it must be related to a breaking of the axial 
symmetry. Let us discuss this using abelian bosonization. 
Here we follow the notation and convention of Ref. ||, 
and take the direction of the magnetization as the quanti- 
zation axis (z-axis.) The breaking of the axial symmetry 
allows a series of operators e 2nmR< ^, where (j> is the dual 
field, R is the compactification radius, and n is an inte- 
ger. If one of these operators appears in the Hamiltonian 
with a non- vanishing coefficient, and it is relevant (in the 
Renormalization-Group sense), we expect an energy gap. 

A simple possibility of axial symmetry breaking is the 
exchange anisotropy, including the dipole-dipole interac- 
tion. However, the magnitude of the anisotropy is of or- 
der of 1% of J [§]|l]]. From a bosonization analysis, which 
we do not discuss in detail in this letter, we estimate the 
gap induced by the exchange anisotropy as 10 -5 J at the 
field of 7 T. This is too small compared to the exper- 
imental value up to 0.3J. Thus we must seek another 
mechanism. 

Cu Benzoate has alternating crystal axes, which gives 
an alternating g-tensor. Due to the alternating (/-tensor, 
a uniform applied field produces an effective staggered 
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field on the spin chain as pointed out in Ref. JjJ. More- 
over, an additional contribution to the effective staggered 
field comes from the DM interaction, which is also present 
due to the alternating crystal axes. Both lead to a trans- 
verse staggered field, namely the direction of the stag- 
gered field is (almost) orthogonal to the direction of the 
magnetization. As we will see, these two contributions 
to the staggered field are of same order and both are 
important in analyzing the angle-dependence of the gap. 

According to Ref. [pi, the local <?-tensor for Cu ions 
is given by g = diag(2.08, 2.05, 2.36) in the local princi- 
pal coordinates. Due to the alternating direction of the 
oxygen octahedra around the Cu ions, the principal axes 
of the g-tensor alternates along the chain. In the ex- 
periment, the field is applied in the principal directions 
(a", b, c") of the total exchange anisotropy. (For details, 
see Ref. §.) 

The <?-tensor in a", b, c"-bases [0 is given by 




2.115 ±0.0190 0.0906 \ 
±0.0190 2.059 ±0.0495 , (1) 
0.0906 ±0.0495 2.316 / 



where ± corresponds to the two inequivalent sites. 
For example, if we apply the magnetic field in in- 
direction, the effective staggered field generated by ex- 
tensor is (0, 0.025, 0)H in a"&c"-coordinates. (The sign 
of the staggered field is defined by referring to the even 
sites.) For field applied in b" and a" directions, it is 
(0.0095, 0, 0.025)7? and (0, 0.0095, 0)H, respectively. 

On the other hand, ignoring other than the nearest- 
neighbor interaction, DM interaction in the chain can be 
written as 

ff DM - 52(-lYD ■ (Sj x S j+1 ). (2) 
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Note that the factor (— l)- 7 is present, as required from the 
crystal structure [ra2]. When a magnetic field is applied, 
an effective staggered field is generated through the DM 
interaction. While it is possible to see this by a Mean- 
Field argument, it can be deduced from the following 
exact mapping. In fact, we can eliminate the DM inter- 
action by a redefinition of the spin variables For 
simplicity, let us assume the D points in the z-direction. 
Then, the Heisenberg Hamiltonian with a DM interaction 
is given by 

H = \ + J*S+Si j+1 + (h.c.)] 

j 

+ J^2i^2j-1^2j + ^2j^2j+l]i (3) 
3 

where S ± = S x ± iS v and J = J + iD. By the rotation 
about z-axis by an alternating angle 

St, - S+e ia / 2 , S+ -> S+e--/ 2 , (4) 



where tana = D/J, the Hamiltonian is transformed to 

3 

+ ^ / Jfffj-lfff? + ^2.7 ^2.7+ ll- (5) 
3 

Namely, the DM interaction is eliminated, resulting in an 
anisotropic exchange coupling. This anisotropy would be 
cancelled by the exchange anisotropy before the redefini- 
tion of eq. (Q), under some assumptions Js) . In any case, 
the resulting anisotropy would be small for Cu Benzoate 
and is neglected in the present letter. 

When an external magnetic field is present, the Zee- 
man term appears in the original Hamiltonian. For ex- 
ample, if we apply the magnetic field in ir-direction, the 
Zeeman term is -Hzeeman = — ^o^-l^J where Ho is the 
external field After the redefinition (Q), the Zeeman term 
gives 

- H cos f^Sy-Ho «n f £(-l)^J. (6) 
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Namely, the effective staggered field of strength 
-ffo sin (a/2) is generated in y direction. For general di- 
rections of the DM vector D and the external field Ho, 
the direction of the effective staggered field is Ho x D. If 
D << J, then a ~ D/J and the staggered field is given 
by H x D/2. 

In Cu Benzoate, a staggered field is already present be- 
fore the redefinition, due to the alternating g-tensor. The 
total effective staggered field is obtained by the redefini- 
tion of the Zeeman term Q together with the alternating 
<?-tensor. For a small alternating part of the g-tensor, the 
total effective staggered field is given by a sum of two ef- 
fects. Both effects produces a transverse staggered field 
(orthogonal to the direction of the applied field), apart 
from the small longitudinal component due to the uni- 
form part of the off-diagonal elements of the g-tensor. 
We neglect the longitudinal component of the staggered 
field, which is actually very small in Cu Benzoate. 

Thus we are led to consider a one dimensional Heisen- 
berg Hamiltonian with mutually perpendicular uniform 
field H and staggered field H 

H = Y,[JSi ■ S i+ i - HSf + h(-iySf], (7) 

i 

with H » h. It is instructive to analyse the model (0) 
in the standard spin-wave theory approximation (lowest 
order 1/s expansion) p], although some of the conclu- 
sions will be modified when we take into account one- 
dimensional quantum fluctuations more accurately. The 
classical groundstate is a canted antiferromagnetic struc- 
ture. The canting angle measured from the x-axis, 8, 
determined by minimizing the classical energy is the so- 
lution of: 2Jss'm28 — Hsin8 — hcos6 = 0. Now con- 
sidering fluctuations around this classical groundstate to 
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lowest order in 1/s gives 2 branches of spin- waves, in the 
antiferromagnetic Brillouin zone (|fc| < tt/2) with ener- 
gies: 



E±(k) = { [2 Js cos 29 + H sin 9 + h cos 9 ± Js(l 
cos fc] 2 - [Js(cos26> + l)cosfc] 2 } 1/2 . 



cos 20) 
(8) 



In the case /i = 0, the minimum energy of two modes 
are E + — H and £L = at k = 0. The - Goldstone 
mode corresponds to a precession of the spins around 
the x axis. A non-zero staggered field h gives this mode 
a finite gap. To leading order in h but all orders in H 
this is given by: 



£L(0) = y/AJsh[l + (H 2 /8J 2 s 2 )][l - (H/4Js 



,2il/4 



(9) 



Note the singular dependence on h but the weak depen- 
dence on H; £L(0) is essentially independent of H until 
H s» O(Js). Conversely, the upper mode E + (0) depends 
only weakly on h but strongly on H. In the case h = 0, 
the existence of the upper mode at energy H is more rig- 
orously established |l(J without the spin-wave approxi- 
mation. The upper mode E+ is presumably observed in 
the experiment: the peak at higher energy Tiuj ~ 0.8meV 
in Fig. 2(a) of Ref. |Q is consistent with the upper mode 
energy H = 7T ~ 0.8meV. Taking into account ID crit- 
ical fluctuations, the power-law behaviour of the lower 
gap B-(0) oc /i 1 / 2 is changed to h 2 ^ 3 , as we will discuss 
below. It is reasonable to expect the weak dependence of 
the lower gap on H to remain true. 

The low-energy behavior of the system should be well 
described by Abelian bosonization. In the bosonization 
approach, the only effects of the uniform field H is shift 
of the Fermi momentum k,F and the renormalization of 
the compactification radius R. The transverse staggered 
field is mapped to the operator cos (2nR<j)). Thus, the 
effective low-energy theory for the model (0) is given by 
the sine-Gordon model with the Lagrangian density 



C = -{d^f + Xcos(2wR(p). 



(10) 



staggered field. If we neglect the effect of the uniform 
field H, d = 1/2 and thus A - h 2 ^ = h om . Precisely 
speaking, there is a log-correction due to the presence 
of the marginal operator [jll): A ~ /i 2 / 3 | log/il 1 / 6 . The 
log-correction is not significant in the present case. While 
the field-theory argument gives the exponent for the gap, 
it does not determine the magnitude. Thus, we studied 
numerically the excitation gap of the S =1/2 Heisenberg 
antiferromagnetic chain under a staggered field (but with 
no uniform field), by Lanczos method up to 22 sites. The 
result is shown in Fig. 1. We found that, for small stag- 
gered field h, the lowest excitation gap to the total mag- 
netization 1 sector behaves consistently with the field- 
theory prediction. We fixed the proportionality constant 
as 



1.85(^) 2 / 3 J|log^r 



(11) 



from the numerical result. 



Smallest Transverse Gap induced by Staggered Field 



gap 0.3 - 




The operator cos (2ttR<P) has dimension irR 2 , and is more 
relevant than the exchange anisotropy. Actually, it is the 
most relevant operator in the system. At zero field, the 
system is isotropic and R = l/y/2n); the dimension of the 
staggered field operator is 1/2. The radius R is affected 
mainly by the uniform field. Its effect may be estimated 
from the Bethe Ansatz solution for the Heisenberg model 
under a uniform field. The dimension of the operator 
cos (2wR4>) is reduced to irR 2 = 0.41 at H ~ 0.52 J (7 
T for Cu Benzoate). While the uniform field does affect 
the low-energy excitation, the effect is not drastic. This 
is consistent with the spin-wave calculation. 

A Renormalization-Group argument shows that the 
gap A is proportional to h 1 ^ 2- ^ where d is the dimen- 
sion of the relevant operator and h is the total effective 



FIG. 1. The lowest excitation gap to S z = 1 sector in 
the Heisenberg antiferromagnetic chain, induced by the stag- 
gered field h. The gap is obtained by an extrapolation of 
finite-size gap by Lanczos method upto 22 sites. Both gap 
and h are measured in unit of the coupling constant J. The 
data is well fit by the field-theory prediction h 2 ^ 3 \ log /i| 1/ ' 6 , 
with a coefficient 1.85. 

Since the effective staggered field h is proportional 
to the applied field Hq, the gap should be propor- 
tional to H 2 / 3 . This is in a good agreement with the 
experiment [Q in which the gap is found to scale as 
i/rj0.65(3) for three directions of the magnetic field. This 
supports our basic claim that the field-induced gap is 
due to the effective staggered field. If we include the 
change in R 2 due to the uniform field, the exponent 
changes to 0.63 at Hq = 7T. Taking an average, the 
agreement with the experiment is improved. Moreover, 
the sine-Gordon model ( |l0| ) is integrable. The elemen- 
tary excitations are given by soliton, anti-soliton and 
soliton-antisolition boundstates ("breathers"), and their 
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exact mass ratios are available Jl2| , |l3[ ] (For an intro- 
duction, see for example Ref. ]14[.) At the isotropic 
point R = \/^/2tt, there are two kinds of breathers, 
and the mass of the lighter breather is degenerate with 
the soliton/anti-soliton, forming a triplet. The mass ra- 
tio of the triplet and the singlet (heavier breather) is 
1 : \/3- When the SU(2) symmetry is broken, the 
triplet is split and the mass ratio for light breather, (anti- 
)soliton and heavy breather is 2 sin [7r 2 i? 2 /(4 — 2irR 2 )] : 
1 : 2sm[vr 2 J R 2 /(2-7ri? 2 )]. For H = IT, irR 2 = 0.41 
and the ratio is 0.79 : 1 : 1.45. This is close to the ratio 
of three peaks observed in Figure 2 of Ref. |], at 0.17 
meV, 0.22 mcV and 0.34 meV. (0.77 : 1 : 1.55). 

In the experiment, the magnitude of the staggered field 
due to the alternating g-tensor depends on the direction 
of the applied field. The ratio is 0.019 : 0.053 : 0.049 for 
field applied in a",b,c" directions. This can be under- 
stood in our theory, because the proportionality constant 
between h and Hq depends on the field direction. If we 
only consider the staggered g-tensor effect, the ratio of 
the gap is 0.019 2/3 : 0.053 2/3 : 0.049 2/3 ~ 1 : 2.0 : 1.9. 
This does not explain the observed ratio of the gap 
0.55 : 1 : 2.0. In particular, the order of gap for b and c" 
is reversed. Thus, it is necessary to include the effective 
staggered field due to the DM interaction, in order to 
explain the gap. 

In general, the magnitude of D is argued to be of order 
of (Ag/g)J where Ag is a shift of g- factor in the crys- 
tal §. In Cu Benzoate, Ag/g ~ 0.1. While more precise 
estimate of D in Cu Benzoate is unknown, it should be 
in ac-plane (or equivalently a"c"-plane) from the crystal 
structure fflM. Thus D is specified by two parameters, 
for example, by D — \D\ and the angle x between D 
and a"-axis. We first determined D so that it reproduces 
the experimentally observed angle dependence of the gap 
a" : b : c" — 0.55 : 1 : 2.0. We found two solutions: 
(X, D) = (0.22, 0.034J) and (-0.0066, 0.10 J). { X is in ra- 
dians.) Both directions are close to a"-axis (or a'-axis) as 
claimed in Ref. Moreover, both values of D are con- 
sistent with the general estimate D ~ (Ag/g) J ~ 0.1J. 
Thus, it has been shown that a reasonable magnitude of 
the DM interaction can give the angle dependence ob- 
served in the experiment. 

We can also estimate the magnitude of the gap from 
our theory, using (p"l|). For the former solution (%, D) = 
(0.22, 0.034 J), the gap for H = 7T applied in in- 
direction is 0.096 J. For the latter solution (%, D) = 
(-0.0066,0.107), the gap for the same field is 0.15 J. 
Both gives a correct order of magnitude compared to the 
experimental value ~ 0.2meV = 0.13 J in Ref. Jj]]. More 
quantitative comparison would require further analysis 
of the specific heat data, since aspects of the treat- 
ment in Ref. Q could be questioned. They fit the low- 
temperature specific heat by six independent massive 
bosons with same gap A, but the sine-Gordon theory pre- 



dicts four elementary excitations with different masses, as 
discussed. While their estimate presumably gives correct 
order of magnitude, the precise value would be changed 
by a refined analysis. Experiments with other field direc- 
tions would provide a further check of our theory. 

Finally, we comment on other consequences of our the- 
ory. From a scaling argument, the staggered magnetiza- 

1/3 

tion behaves as H . The direction of the staggered 
magnetization is given by the effective staggered field. 
Thus our theory could be tested if the staggered magne- 
tization is measured. Moreover, by the redefinition (Q), 
the physical spin operator corresponds to a rotated spin 
operator in the Heisenberg antiferromagnet without the 
DM interaction. While it has no drastic effect on the 
neutron scattering experiment, it affects the susceptibil- 
ity measurement dramatically. The observed susceptibil- 
ity Xcxp is given by a linear combination of the uniform 
susceptibility %„ and the staggered one \s of the Heisen- 
berg model. Since the latter diverges at low temperature, 
Xcxp would also diverge at low temperature. This could 
explain the enhancement of the susceptibility observed 
in Ref. p"5| , though a quantitative theory would require 
inclusion of interchain interactions. Further discussions, 
including details of the arguments in the present letter, 
will be given in a future publication. 
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